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Q ; 1 Introduction 

\0 ■ 

f^) ' Granular materials partake almost dramatically at times of the properties of 

C*~) , solids and, under different circumstances, of some properties of gases. Some 

scientists have suggested that they represent a new state of matter, to be treated 
. per se; because of their ambiguity the search of an overall continuum model for 

Oh' these media is quite intriguing. Even if one confines inquiry to fast flows of 

fH . sparse granules (as in granular gases, see e.g. [2]), the range of phenomena that 

+^ ' may occur and that should be described by an adequate theory is vast. Within 

S this area one is lured into adopting such appealing terms as gross granular heat 
and gross granular temperature, though the appropriateness of the borrowing 
. can be challenged as the definitions of those gross quantities involve kinetic 

entities, whose detailed evolution can be in part ascertained rather than being 
totally unpredictable on principle and, at most, appraised through statistics. In 
any case, when exploiting vague analogies with the thermodynamics of perfect 
gases, one must make allowance for the fact that, in generic granular flows, 
speeds of agitation need not be, as they are in gases, many orders of magnitude 
greater than the average velocity or the velocity of any reasonably behaved 
observer. Thus one must be wary of a naive transfer to each moving element 
in a continuum of properties inferred from experiments within containers fixed 
to the walls of a laboratory; questions of objectivity become important if not 
paramount. One way to avoid difficulties would be to evaluate the relevant 
variables from some intrinsic local frame. The traditional way in continuum 
mechanics is to appeal to gross position and velocity gradients, a way that 
cannot be satisfactorily pursued in our context. Notice that akin problems 

"The paper completes, extends, hopefully makes more perspicuous topics covered by an 
earlier report of limited circulation [1]; the report itself was an excerpt of a talk given in 
Bologna in 1999. 
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arise in the analysis of vibrations of large and very flexible space structures (or, 
more down to earth, of the motion of a misshaped soap bubble in wind): the 
centre of gravity need not be superposed to a place occupied by some element 
of the structure; even if perchance such superposition occurs, the movement in 
the vicinity of the centre might have little to do with the gross motion of the 
structure; hence the need of an alternative notion of gross rotation. A solution 
is easy to find: the set of balance equations of momentum, moment of inertia, 
moment of momentum can be interpreted as the tools that allow us to figure 
out the background against which an objective measure of residual energy of 
agitation can be achieved, with consequent satisfactory portrayal of 'thermal' 
concepts. 

Again, the many different scales of events possible in a generic granular flow 
entail a reconsideration of the ties of the concept of temperature with the ki- 
netic circumstances of ensembles of granules having diverse (and not necessarily 
canonical) distributions of energy, thus including easily, for instance, negative 
and tensorial absolute temperatures. 

Finally, when one pursues a continuum model of granular flows, one must 
abandon the fundamental tenet that it be possible to define, that is to iden- 
tify for ever, each material clement. Then test domains of minute diameter 
and the corresponding local averages take a central role; the question how to 
average near the frontier of the body and thus how to model boundary effects 
must be addressed. Moreover, in the formulation of constitutive relations, valid 
cither in the interior or at the boundary of the body, local space and time- 
correlations might take precedence over material space-gradients and material 
time-derivatives. 

However, whichever changes be necessary, ultimately some alternative paradigm 
should emerge and it is attractive to note, from a cursory review of some pro- 
posals for extended versions of continuum mechanics (e.g., the theory of hy- 
poelasticity) or continuum thermodynamics, some stunning similarities in the 
evolution equations which are arrived at. True, those theories all presume the 
absence of local mesospin (which may be instead relevant in granular flows). 
However, next to the equation of continuity and Cauchy's equation an addi- 
tional equation concerns the shuffling motions; it involves a stirring tensor and 
rules the evolution of a Reynolds' tensor, a symmetric tensor which both in hy- 
poelasticity and in extended thermodynamics coincides with momentum flux. 

Precisely, the theory of hypoelasticity [3], Grad's theory of 13 moments [4], 
extended thermodynamics [5], Jenkin's theory of fast flow of granular materi- 
als [6], all suggest the addition to the classical equation of balance of another 
equation which rules the evolution of the symmetric Reynolds' tensor H: 

A 

pH = Si — div s + Se', 
at the same time, for Cauchy's stress tensor T, a very simple constitutive law 

A 

is suggested: T = pH . p is density, H the convectcd time derivative of H, s 
is the third order stirring tensor; the tensor Si represents internal equilibrated 
actions, and Se the external stir. 
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Here, within the mechanics of mass points, an elementary analysis which 
involves predominantly velocities rather than places, is shown to lead to a global 
equation of motion which suggests, by analogy, within continuum mechanics, 
just the last balance equation mentioned above. Precisely that equation may 
be relevant, in general, when modelling kinetic bodies (i.e., bodies which are in 
a permanent state of flux and possess no natural paragon placement) and, in 
particular, granular gases. 

2 Adscititious topics in elementary mechanics 

Consider a system S of N mass points; let x (r) the current centre of gravity 
of S and /i the total mass; //W, the mass of the i—th point; a;W ( T ), its current 
place; y^ , its position vector with respect to the centre of gravity in a paragon 
setting (e.g., the initial placement). Take for R(t) a proper orthogonal tensor; 
split x^ into the following sum: 

x {i) (t) = x (r) + R (t) s (i) ; «W = jfW + ( T ) (i) 

involving a global rigid displacement (x(r),R(r)) and an individual shuffle 
vector 

In a motion the velocity of each point is given by the sum 

±M ( T ) = x (r) + e (q (r) ® y« (r)) + (r) s« (r) (2) 

of speeds of entrainment (the first two addenda) and of agitation (the last ad- 
dendum): here — x^ — x, e is Ricci's permutation tensor, q is the speed of 
rotation so that RR^ 1 = eq. 

Differentiating again (J2J) with respect to time one can exhibit explicitly en- 
trainment, relative, Coriolis' components of acceleration: 

i w = x + e (q $ y«) +? 2 (J-c®c) y w + i?s w + 2e (g ® i?s w ) ; (3) 

here c is the unit vector associated with q, i.e., q — \q\ c, and / is the identity 
tensor. 

The total kinetic energy /in is given by 

= + ^w 2 (/ — c ® c) • y + 

5 E ^ } 2 + E ^ w (> w ) • e (« ® . w 

where /iY is Euler's inertia tensor 

A1 y = ^ M «y«®y«. (5) 
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Notice that \iJ ', J = (trY) I — Y , is the usual tensor of inertia, and that 
n(I — c® c)-Y = c-Jc coincides with the moment of inertia around a baricentric 
axis parallel to c. 

Actually, one could have started this elementary investigation taking for x 
any moving point and enquiring subsequently how x (r) should be chosen so that 
t^x 2 differ the least from the total kinetic energy; the answer would have been 
to choose x so that it coincide always with the speed of the centre of gravity. 
Similarly, one can choose q in such a way that the speeds e (q (g) j/W) best fit 
the motion relative to x (r), in the sense that the discrepancy in kinetic energy 
is minimal 

1 



2Z ^ ~(q®y (l) 

i 

one needs only take for q the solution of the equation 



(6) 



Jq = k, (7) 
where fik is the moment of momentum with respect to the centre of gravity 

M* = J> (<) y w xyW. (8) 

i 

With that choice of q the last term in J3J cancels out. The kinetic energy is split 
into the sum of the two, observer-dependent, usual terms as in a rigid motion 
and an observer-independent term due to agitation. 

To prove the statements above notice that the last term in Q could be 
written successively as follows: 



yH) X 



y W x (g x y {1 ' 



whereas the second addendum in the right-hand side of J1J is exactly the con- 
tribution to fiK of a rigid rotation around x. 

Knowledge of x and q (and hence of a; and R) offers the chance to fix a 
background against which the motion of agitation can be measured objectively. 
To determine x and q, appeal can be made to the global equations 



fix = f 



(9) 



where / and m are resultant and resultant moment of external forces on S. 
Success is obvious if there is no agitation and the motion is trivially rigid. 
However, in general one is not so fortunate because there may be an influence 
of shuffle and agitation certainly on J and, perchance, on / and m; notice, in 
particular, that J depends critically on expansion/contraction of S because 



J = trY) I -Y 



(10) 
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and 

Y = 2symK, (11) 

where 

^=xyv i) ®2/ w (12) 

is the tensor moment of momentum, the evolution of which must be known to 
evaluate the evolution of Y. In a rigid motion K depends on Y and q only, 
K = — (eq) Y , and no further developments are required. Instead, in general, 
as remarked above, a deeper analysis of the motions of S needs be effected in 
advance. We need to modify the developments at the beginning of this section 
as a premise; rather than Q take now 

x ® (r) = x (r) + G (r) s {i) (r) , (13) 

where G is a tensor with positive determinant; correspondingly, 

iW = x + By® + Gs {i) , B = GG~ 1 , (14) 

and 



W =£ + (b + B 2 ) y^ + 2BGs^ + Gs {l \ (15) 



if is equal to Gs' 1 '. 

Thus, shuffle and agitation are now supposed to be remainders beyond an 
affine, rather than rigid, motion. The kinetic energy is given by the sum 

\^ + \p*r (BYB T ) + \ ]T M W (Gs«) * + ]T ^ (g*») ■ (b V ®) . (16) 

i i 

However, if B is chosen so as to make /iK — fj,^ (i??/ 2 )) 2 a minimum, i.e., as 

i 

to satisfy the condition 

BY = K T , (17) 
the last term in (|l(j|l vanishes; in fact, even the mixed kinetic tensor 

2^WflW®ByW (is) 

vanishes, not only its trace, so that a relatively compact expression is available 
also for the kinetic energy tensor [iW = | follows: 

M W = -fix <8> i + \^iBYB T + i/xff. (19) 

where, if H* is Reynolds's kinetic tensor of agitation, H is its transformed 
expression by G: 

/iff* = ® . # = Gff*G T . (20) 
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It is important to notice, in particular for later developments, that all vectors 
s^' and the tensor are absolute, i.e., not affected by changes or movements 
of the observer. 

To determine x and B, appeal is now made to the first one of © and to the 
equation of balance of the tensor moment of momentum 

fj, (K -BK-H) = M - A, (21) 

where M and A are respectively the tensor moment of external and internal 
forces acting on S 

M = y (l) ® f exm ' A = -J2y {z) ®/ mt(l) - (22) 

i i 

Because the vector moment of internal forces vanishes, A is a symmetric tensor. 
Remark that can be read as the assertion that the 'reference' tensor of 
inertia Y* = G~ 1 YG~ T is constant irrespective of shuffle. 

We still need the evolution equation for H; it is easily obtained from New- 
ton's law, multiplying both members tensorially by , summing and trans- 
forming through G: 

fi (H + BH + HB T ^j =S-Z; (23) 
here S is the stirring tensor of external forces 

S = 2symJ2 ( G ^ W ) ® (24) 

i 

and a similar definition but with opposite sign applies to Z , involving internal 
forces. The queer choice of sign for A and Z has to do with a convention 
appropriate in a distinct, later context. Finally, if one remarks that 

1 1 
K = symK + -ere , M = symM + 2 em ' (25) 

the system of equations to explore becomes 

= /, 
k = m, 

( fisym (K - BK} = symM + fj,H - A, ( 26 ) 
fx (h + BH + H B T ) = S - Z. 

Strictly, the first two equations bear upon the preferred reference, if only its 
rotational speed q is defined through the relation 

q = ,r x k, 

whereas the last two lead to global hints on agitation; but, generally, the two 
sets are strongly linked. Actually, for some later purposes, it is more convenient 
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to keep together the two equations for the symmetric and skew components of 
K and to add to the list the evolution equation for Y: 



V>x = f, 

n(k- BK \ = M - A + fiH, 

Y = YB T + BY, 
fx (h + BR + H B T ) = S - Z. 



(27) 



Notice, in the second and fourth equation, that the quantities between brack- 

A A 

ets express the convected time derivative K of K and H of H, based on the 
spin tensor B. 

If all sW vanish, then all sW maintain their initial values, H vanishes and 
the local motion is affine, or pseudo-rigid in the therminology of Cohen and 
Muncaster [7]. Their results could be borrowed here; one needs to study only 
the reduced system of the first three equations of (|27Jl in x and B. 

More generally, if /, M, A, Z and S depend at most on x, x, G, B and H, 
then Il2(ill can be interpreted as a differential system in x, G, and H. 

Circumstances could be called kinetic when the system has no physically 
relevant paragon setting, its behaviour is ruled by abrupt responses to current 
circumstances and, as a consequence, in our model, /, M, and S need depend at 
most on v = x, B and H ; then (|27[1 becomes a first order system in v, B, Y and 
H. This system would merit scrutiny on its own, perhaps under appropriate, 
special choices of /, M, A, Z and S. 

A case of special interest is met when A = fiH or, at least the quantity 
(A — fiH) together with / and M depend at most on v, B and Y; then the first 
three equations can be dealt separately from the last. 

More subtle is the case in when /, M and S depend also on k+^ ; where cpj is 
the potential energy of internal forces, supposing that they all be conservative; 
or, rather, when they depend also on the manner the mass points of S can 
be classed in families of increasing total energy, or, more deeply, in families of 
approximately equal Reynolds tensor. 



3 The kinetic energy theorem and some corol- 
laries 

A tensor kinetic energy theorem can be derived easily from H27JI by adding term 
by term first, second and fourth equation after tensorial multiplication of the 
first by x, of the second by B, and after multiplication by | of the fourth; finally, 
by taking the symmetric parts of all terms: 

fxW = sym [ix ® x + fi [bYB t + B 2 YB T ^j + ^ f iH = 

= sym [fix®f + B (M -A)] + ~(S-Z). (28) 
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More particularly, but also with deeper meaning, if one takes the trace, one 
arrives at the more usual kinetic energy theorem 

k = x ■ f + (M - A) ■ B T + -tr (S-Z). (29) 

Remark that 

M-B T + \trS = Y J V i -f i E, 
-{A-BT + \trZ)=^-n- m 

i 

The standard requirement that the power of internal actions be invariant for 
any rigid change of speed, when B is an arbitrary skew tensor, also leads to 
the condition 

A G Sym, 

which was already noticed on equivalent grounds. 

If potentials tp E , tp I exist for external and internal forces, respectively, then 

fh = ^ ' /1 = t?"' (31) 
a theorem of energy conservation follows: 

K + (p = const , ip = (p E + (pj. (32) 

Notice that, on the one hand, for Q14[). 

dip 

tp- x 

and, on the other hand, 



E5-(*+^^ (0 ). 



dx ^ &cW ' dG ftcW ' 

so that 

d< P \ , n I dl P « „fOl ^ dip 



Juxstaposition with the right-hand side of (|29|l suggests the 'constitutive laws' 



0a« j ' 



or, more precisely, 



S = 2sym s W 



d<PE 
<9s« 



tp E may be influenced also by x and skwG, besides syraG and sfi'\ on the 
contrary, tfj must be observer independent, hence it must not involve x and 
may depend on G only through the product C = G T G, so that 



(34) 



Of course, in general, A and Z have also dissipative components beside the 
conservative components expressed, as above, through ip T . 

Sometimes there is an interest for a 'reduced' theorem involving the 'gross' 
kinetic energy tensor W per unit mass 

1 1 rp 

W = -x®x + -BYB T . 
The theorem is signified by the reduced equation 

fj,W = S ym {fix® f + B (M - A)} , 

and leads, by difference from l|28|l (and when the trace is taken), to a 'principle' 
of energy balance. The latter is, in the present context, nothing else but a 
corollary of the last equation H27|l but with a different promotion of terms, the 
leading role being played by the 'internal energy' 



e = -trH 
2 



ipj. 



(35) 



As we can avail ourselves of the more powerful relation J3SJ|, we do not 
pursue here the consequences of that corollary. Actually, for later purposes, one 
can write the last equation l|27|) in an equivalent form which approaches the 
principle of conservation more closely, using a tensor of energy 



E=-H + -ififl: 
2 V 1 ' 



precisely 



where 



n (e + be + eb 



f 2 
Z - ~<p T I - -ifi I symB 



-2sym 
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4 Energy distributions 

The topics of this section are textbook affairs; because of the interest here in 
non-canonical instances, they are recalled nonetheless in essence and with the 
appropriate slant. 

4.1 The scalar case 

When the mass points are very numerous, though with bounded total mass //, 
the 'averages' x, k, Y, etc. acquire prominent import. At the same time, as men- 
tioned at the end of the previous section, the resultant actions on <S, expressed 
by /, M, S, may come to depend (not only on kinematic variables such as v, 
B, Y, H, but, as hinted, also) on the way mass-points can be parcelled out in 
families, each family comprising points with an energy of agitation falling within 
a limited range, say [(j — 1) /isd, j /j,eS], j = 1, 2... (here S is a positive number). 
The fraction ^ of mass points belonging to the j— th family is measured by a 
non- negative constant 7-. 

Then an histogram can be drawn as a graph of a piecewise constant function 
7 (£) having the value 7 • for £ within the interval [(j — 1)6, jS] . Notice that there 
will always be a value, say J of the index, such that Nj > 0, whereas all Nj 
with index larger than J vanish; that it is so because the total energy me is 
(approximately equal to and) not less than fieS (j — 1) Nj-, that each term in 



the sum cannot exceed cP ; that any non-null value of Nj is a positive integer, 
hence no less than 1; and, in conclusion, that J — 1 < 6~ . 
The function 7 (£) satisfies a normalization condition 




J o 7(o« = *ETi-i = ]vE^ = 1 



3 3 



A second important relation is derived easily; remark that 




and that, on the other hand, 




3 



3 



or 



1 <<5 2 ]Tj7, <! + *• 



3 



In conclusion. 
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Our analysis becomes more fluent if we proceed to smooth out the histogram 
(perhaps imagining that the range of 5 are taken ever smaller) to become the 
graph of a (continuous, even smooth) function with the properties that 

7 (£) d£ gives the fraction of mass points with energy of agitation within the 
interval (/i££, jJ,e (£ + c££)) and that the following normalization conditions apply 

/•oo />oo 

/ 7(0^ = 1 , / ^7(0^ = 1- 
Jo Jo 

It is easy to contrive distribution functions satisfying all conditions noticed 
so far for 7. Take any function A (£) defined over [0, +00) with non-negative, not 
everywhere null values (even a measure) and integrable, together with £A (£) , 
over [0, +00) 

/>oo />oo 

/ \(l)dl = p , / l\(l)dl = a; 0</>,<7<oo. 
Jo Jo 

By choosing 

one obtains just one of the desired functions. 

It is an easy matter to check that the following choices for 7 satisfy all 
requirements mentioned above. It is appropriate to emphasize that the abscissa 
£ for the histogram is chosen here as to be non-dimensional and such that 
mass points for which £ = 1 have energy per unit mass exactly equal to the 
total energy per unit mass. The use of non-dimensional variables may give an 
impression of excessive specialization; in the formulae, in fact, the choice of 
constants is mandatory: 

(i) Canonical: 

7 = e" £ , for all £. (36) 

For some systems, such as monoatomic gases, this distribution is requisite 
under 'quasi-static' conditions. 

(ii) Power law: 

7 = 24(2 + £)~ 4 , for all £. (37) 

Remark that other negative powers different from —4, strictly less than 
—2 could be contemplated; the numerical factors would then be obviously 
different. 

(iii) Piece- wise constant: For any constant (3 G [0, 1), take 

7 (£) = 0, when < £ < (3 or £ > 2 - /3, 

7(0 = 2(1-/?) ' Wh6n P<t <2 ~P- ( 38 ) 
The limit for (3 —> 1 is a measure with an atom at £ = 1. 
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Piece- wise linear: For any f3 £ [|,3], take 

7 (O = 0, when £ > /3, 



7 (£) = 2/T 3 [3 (2 - /3) £ + (2/3 - 3) /?] , when < £ < /?• (39) 

Piece- wise exponential: Given the constants /3, positive; £ x > 0, £ 2 > ^; 
and a real; take 

7 (O=0, when < £ < £ x or £ > £ 2 , 

7(0=^"^, when SGKi,^], 



where 



and a is a solution of the equation 



(1 + a£ 2 ) e~ a ^ 2 - (1 + <i) e~" €l 



(40) 



To lighten developments, we refer below only to the choice O = 0, when 



7(0 = ^-^, (41) 



where a satisfies the equation 

"0 



1 - e -«« 2 



(42) 



An inspection of l|42|) requiring only accurate evaluations of orders of mag- 
nitude shows that: There is one and only one value of a satisfying l|42|l 
for each choice of £ 2 larger than 1. There are no values of a satisfying 
(|42|) for < 1- The function a(£ 2 )j thus defined, is strictly increasing 
from — oo to 1; it vanishes for £ 2 = 2; its approximate expression in the 
neighborhood of 2 is 

oc * \ - 2) . 

When £ 2 tends to 1 and a to — oo, the distribution (|41|l approaches a <5— 
function with an atom at £ 2 = 1. When £ 2 tends to the value 2 and a 
to zero, the distribution l|4*T|) tends to be piece-wise constant: 7 (£) = \ 
for < £ < 2 and null otherwise. When £ 2 tends to 00 and a to 1, the 
distribution tends to be canonical. This distribution is reminiscent of one 
which is suitable for quantum systems though allowing in that case only 
a finite number of states. As is well known, for them, with appropriate 
care, a sudden transit from a distribution with a positive into one with a 
negative can be achieved experimentally. 
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(vi) Sinusoidal: Given a constant a less than ^- take 

4a 



4 

Si , \ ' l 

7 (£) = 0, when £ > 2 1 



7T 2 



7({ ,_l( 1 _5)( 1 + aeo .l( 1 .5) { ). (43) 

The limit when a goes to zero coincides with the limit case (iii) when 
[3 — > 0. 7 either decreases or increases with increasing £, depending on the 
sign of a. This distribution may be of interest when energy and numcrosity 
are functions of an angle from a given direction. 

(vii) Fermi and Bose. For j3 non-null 

7 (f)= e l ' 3| «(e' 3 -l)~ 1 -l ~* (44) 



satisfies the first normalization condition; now choose (3 so as to fulfill also 
the second one, which requires that 

1 = Dilog(e' 3 - 1) , (45) 

where 



k=l ■ " 



There are two solutions, one with (3 negative (Bose-Einstein subcase) and 
the other with j3 positive (Fermi-Dirac subcase): 

(3 = -0.814651... , (3 = 1.405050... (47) 



4.2 The tensor case 

As hinted repeatedly, a classification of mass points in families with specific 
energy is inadequate at times to cast satisfactory constitutive laws completing 
system <|26|) or 1)2 7}. One can explore broader alternatives (see, e.g., [8]), in 
particular one can attempt the classification in tensorial terms; i.e., when ip T 
vanishes, in terms of the Reynolds tensor. Precisely, to pursue here the latter 
(special, but most important) case, think of the linear space of symmetric tensors 
and, within it, the manifold N of all rank-1 tensors, i.e., of all tensors N of 
the type v (£> v with v any vector. With 7 a function of the tensor variable 
H~zNH~2, call now 

7 (H~^NH~^ (H- 1 ■ dN) 

the number fraction of mass points with Reynolds tensor in the immediate 
neighborhood of N; imagine again the histogram of 7 smoothed out to render it 
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continuous and differentiable. It would also be integrable over M together with 
its product by H~ l N and with the properties of normality 

J "f(H-^NH-^H~ 1 -dN = l, 

J (H- 1 ■ AT) 7 (h-^NH-^ (H- 1 ■ dN) = 1. 

Then, formally, one can proceed in analogy with the scalar case; in particular, 
one can introduce the extended canonical distribution, valid on A/", 

7 = e 

(see, e.g., [9]) and explore also alternative distributions. 



5 Granular temperature 

When the intimation at the end of Section 3 is mandatory, knowledge of the 
distribution 7 becomes essential; then the need would appear, from the develop- 
ments of Section 4, to seek another equation (a sort of Boltzmann equation) to 
describe the evolution of 7. Actually it often occurs that either the class of the 
accessible distributions is restricted on physical grounds and each distribution 
is then identified by few parameters, or, at least, only a few variables linked to 
the distribution are essential and suffice. Such is the granular temperature i9, 
which is defined below: 

It is a constant parameter in case (i). If the case (iii) is taken as the limit case 
(v) for a — > 0, the associated temperature is infinite. In case (v), subcase £1 = Q, 
then 1? = a -1 ; hence the graph of temperature against £ 2 shows 1? decreasing 
from to —00 as £ 2 ranges from 1 to 2, and then decreasing from +00 to 1 as £ 2 
grows to +00. It is noteworthy to remark that null temperature can be achieved 
only through negative values; transition from positive to negative temperature 
can occur only through 00, as is known experimentally. 

In distributions different from exponential, as in case (ii), the derivative l|48|l 
is not a constant; one can declare that temperature is meaningless for such cases 
or be satisfied with an average value d w of the derivative over its support. 

More precisely, suppose that, as required for invertibility, 7 (£) be strictly 
monotone and positive only over an interval (£1)^2) e then for the 

weaker definition § w of the temperature, one obtains 
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Summing up, for the special distributions above, the temperature turns out 
to be 

(i) 0=1 ; (ii) <?„, = - ; (iii) = oo ; 



2 

^ *~ = \j=2> H -4a) a ' 



1 



In case (v), subcase £ x = 0, i? = a -1 , hence i? varies with the choice of £ 2 as 
was described above. 

As hinted at the end of Section 2, the temperature (or its reciprocal, the 
temperance), even when allowed to roam over the entire real axis, may be yet 
inadequate to allow closure of system i|26[l , l|27|l for a sufficiently wide class of 
bodies. Then one can explore broader alternatives (see, e.g., [8]) and insinuate 
the idea that a tensor related to the distribution of shufflings may be of help. 
Precisely, one may think to the tensor case discussed in Section 4. Then, over the 
support of 7 on J\f also log 7 becomes meaningful and the following temperature 
tensor can be defined: 

(50) 

As already mentioned, by far-fetched analogy with standard cases, one may 
dream up occurrences when H , S, and Z are functions of the temperature tensor 
on physical grounds; then, the last equation (|27(l becomes an evolution equation 
for the temperature tensor. Naturally such dreams need by substantiated by 
offering some concrete cases where the fantasy bears fruits; for the moment, 
let rest us content to have drawn attention to a possible avenue which is wider 
than that opened up by the mere scalar temperature and hence perhaps more 
attuned to the study of granular materials. 



6 Transition to a continuum model 

The system (I27|l is already written in such a way as to suggest the possible 
evolution equations for a continuum, where the image of the body is a fit region 
rather than a discrete set in Euclidean space. 

Now mass density p takes the defining role for mass; p and the current place 
x, the current gross local shape G, the current tensor of moment of inertia Y, 
and the kinetic tensor H, are all fields. 

The conservation equation for mass is the classical one 

p + pdivx = 0, (51) 
and its twin for tensor moment of inertia is still the third of Il27() : 

Y = YB T + BY , B = GG~ l . (52) 
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The evolution equations for x, G and H can be formally adapted from the 
first, second and fourth of l|27|) : 

px = f, 
p (j( — BK - H^j = M — A 

p(H + BH+ HB T ^j 

The sources per unit volume of momentum, moment of momentum and stir 
need be still specified. The choice is most difficult because those sources might 
not be necessarily of purely local character; at least weakly non-local effects may 
be expected. Thus a vast field of inquiry is opened, possibly, even necessarily, 
involving concepts recalled in Section 5; below reference is made only to very 
special choices advanced in important contexts. They all presume that the gross 
local shape and orientation be invariant: B = 0; Y, constant. This presumption 
implies the irrelevance of equation (|52[1 and the dropping of the second equation 
in system ()41|l. thus implying that A = pH; besides, for /, the standard Cauchy 
form is postulated 

/ = pb + divYi, 

requiring a volume action b and a Cauchy stress S; finally, for £, the special 

constitutive relation £ = pH is accepted. The choices of S and Z remain open 
and again an extended Cauchy type Ansatz is called upon with a volume stir 
and a shuffle flux expressed with the help of a hyperstress; for both stir and 
hyperstress appropriate constitutive relations are finally chosen. 

All these special choices seem to be satisfactory within the particular con- 
texts where they are advanced. But perhaps some general rules for /, M , A, S 
and Z should rather be sought and corresponding fairly general subchapters of 
continuum mechanics exhaustively written. 
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